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A Comment [l[ to our recent article 0] , expressed crit- 
icism on our estimate of the maximum magnetic field 
that can exist inside a neutron star and on the pressure 
anisotropy we found for a magnetized gas of fermions. 
These two points are relevant for the physics of neutron 
stars. With our reply we attempt not only to address 
the Comment's criticisms, but for the sake of the reader's 
understanding, we also try to clarify the connection be- 
tween our and some literature's results mentioned in the 
comment. 



I. MAXIMUM FIELD STRENGTH 



this simple example was to illustrate that the maximum 
field estimate can change just by relaxing the assumption 
of uniform mass and field distributions that was used in 
0. In the second paragraph following Eq.(lO) we wrote 
"we make no claim that this ad hoc model correctly de- 
scribes the real way the field varies with the radius in a 
hybrid star" . More importantly, this ad hoc model was 
not used for any of the main derivations of the paper, so 
it has no relevance for any of the paper's main results. 
Notice that there are other attempts in the literature Q 
of parameterizing the variation of the magnetic field with 
the stellar radius. 



As known, the inner magnetic fields of neutron stars 
are not directly accessible to observation, so their possi- 
ble values can only be estimated with the help of heuristic 
methods. A widely accepted estimate of the maximum 
star's magnetic field ~ lO^^G was done in Ref. us- 
ing a macroscopic analysis based on the equipartition 
of the magnetic and gravitational energies (Virial The- 
orem), and assuming a constant density and a uniform 
magnetic field throughout the entire star. As stated in 
the Comment [l[ , this estimate is in agreement with nu- 
merical simulations that used various EoS of nuclear mat- 
ter 0,0, although it should be mentioned that the results 
of Ref. U were later criticized @ because in their cal- 
culations the authors omitted the compositional changes 
in the EoS due to Landau quantization. Same criticism 
is applicable to Ref. since the magnetic field was not 
considered in the EoS they used. Ignoring the field in the 
EoS of nuclear matter is not consistent at strong fields 
because the EoS is substantially softened by the magnetic 
field for fields well below the proton critical field ~ 10^° 

In Ref. [2] we argued that if the star is a hybrid star, 
meaning its core is made of deconfined quark matter, 
while its outer layers are made of nuclear matter, the pre- 
vious inner field estimates based on a constant mass dis- 
tribution may not be reliable because in the hybrid star 
a much larger difference should exist between the matter 
densities in the outer layers and in the core. In Section 
IIB of our paper, we did two independent estimates of 
the maximum magnetic field of a hybrid star assuming 
a very dense core. The Comment criticism pertains to 
the first of these estimates. Like in [sj, our first estimate 
was based on the equipartition of the gravitational and 
magnetic energies in a star of given mass and radius. The 
main difference with respect to Q was that we consid- 
ered an ad-hoc model with nonuniform matter and field 
magnitude distributions. The only purpose we had with 



The Comment failed to mention that in [2[ we did a 
second, completely independent, microscopic estimate of 
the maximum inner field in the last part of Section IIB 
that not only is model-independent, but it is based on 
the very natural assumption that the magnetic energy 
density should be at most as large as the baryon energy 
density. What this second estimate indicates is that if 
the baryon density in the core is large enough to dccon- 
fine the quarks (chemical potential fi > AOOMeV), the 
maximum magnetic field there can be as large as lO^^G. 
Interestingly, this second estimate is similar to the one 
considered in Ref. 14] , which is one of the papers cited in 
the Comment. In the criterium to estimate the max- 
imum field at the center was based on the stability con- 
dition given by the cancelation between the nuclear fiuid 
pressure and the magnetic pressure, i^matter = B^/2,{see 
Sec. 4.2 in Ref. equivalent to the condition p|| = 
of our paper. Now, taking into account that the main 
contribution to the fiuid pressure is ^ and that the 
chemical potential in the case of quark matter can be 

~ 400 MeV, one straightforwardly obtains B 10^" G 
for the case of quark matter. The smaller field estimate 
^ lO^^G of [3, @1 is connected to the fact that they con- 
sidered nuclear matter instead of quark matter. An ad- 
ditional confirmation that the condition ilmatter = /2 
is reached at field values of higher order than lO-'^^G for 
quark matter, was obtained by numerical calculations in 
Ref. Q, where the EoS of the color superconducting 
phase of quark matter was investigated. 



In summary, the key issue to estimate the maximum 
inner field is how large the baryon density at the core of a 
neutron star is, because its magnitude will determine the 
value of the matter energy density, which in turn sets the 
upper limit for the magnetic energy density achievable at 
the core. 
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II. PRESSURE ANISOTROPY 

In the second section of the Comment [l[ , the authors 
criticize the anisotropy that arises in the calculation of 
the parallel and transverse pressures in the presence of 
a strong magnetic field. According to them, the finding 
of an anisotropy in the pressure and hence in the EoS of 
the magnetized system 0, |^ is inconsistent, because it 
ignores the effect of the Lorentz force density related to 
the magnetization currents. To sustain their claim, the 
Comment authors cite Ref. and also give two ad hoc 
examples. 

The arguments used in the Appendix of [loj and then 
repeated, without any new theoretical development, in 
Section II of sustained that the pressure splitting 
should be eliminated by the creation of an extra term 
related to the work done by compressing the fermion 
gas in the transverse direction against the Lorentz force. 
The new term involves the magnetization current density 
(V X M) X i?, which would make a contribution, MB, 
to the transverse pressure, hence regaining the pressure 
isotropy. Neither the Comment, nor Ref. [l^ offered 
a rigorous proof of such a claim, as the exact Lorentz 
force required for the cancelation was not derived in a 
consistent way, but introduced by hand. 

We highlight that a magnetization current density 
~ V X M is absent in our paper 0, as well as in Q, 
because V x 7\/ = for the static system with uniform 
magnetic field we studied. Assuming a uniform magnetic 
field in our microscopic formulation is well justified, be- 
cause despite the magnetic field strength of the neutron 
star may well vary from the core to the surface, the scale 
length of the field variation in the star is much larger 
than the microscopic magnetic length Q. 

The authors of [l| also tried to sustain their claim by 
considering two ad hoc examples that have nothing to do 
with the problem studied in our paper. The first example 
consists of a system with boundaries. It not only does 
not apply to our study of the EoS of the bulk material, 
but the analysis is done assuming exactly the hypothesis 
it intends to prove, i.e. that the two "total" pressures 
should be equal to guarantee the isotropy in the pres- 
sures. In the second example, they propose to introduce 
in the calculation of the thermodynamic potential, the 
gravitational acceleration, arbitrarily chosen to depend 



on z, in order to produce an inhomogeneous matter den- 
sity distribution and hence a nonzero V x M, because 
M ~ would depend on z. This is a completely ar- 

tificial example, but even worse, it introduces the effects 
of gravity in £7, when it is well known that they are in- 
significant at the microscopic level and hence should be 
ignored in the EoS. Therefore, none of these examples 
provides any reliable proof of the Comment's claims. 

We point out that the finding of the EoS is based on 
microscopic calculations that pertain to the properties 
of the bulk. Energy and pressure can be found in this 
context as the statistical average of the components of 
the energy-momentum tensor. The essence of the origin 
of the pressure anisotropy is contained in the Eq. (81) 
of our paper Q. In this equation we wrote the average 
of the energy-momentum tensor in terms of all the avail- 
able covariant structures of the system. For example, the 
breaking of the Lorentz symmetry by the density gives 
rise to the well-known contribution to the energy, /lA^, 
coming from the second term of (81). In a similar fash- 
ion, the presence of the magnetic field, which breaks the 
spatial rotational symmetry, yields a new tensorial struc- 
ture reflected in the last term of (81). This last term is 
responsible for the pressure anisotropy. 

To conclude, we underline that Ref. [l| has not pro- 
vided any valid proof that the finding of the pressure's 
anisotropy of our paper is inconsistent. Our result was 
rigourously proved in Q. Whether the anisotropy can 
or not be compensated by the currents generated by the 
rotation of the star, a surface magnetization produced 
by boundary effects, etc, is still unclear, but unrelated to 
the microscopic calculations performed in 0] . To answer 
this question, a rigorous macroscopic analysis needs to 
be done, one that uses the equations of magnetohydro- 
dynamics and that is based on a self-consistent treatment 
of the equations of the stellar structure in a strong mag- 
netic field, i.e. using axisymmetric TOV equations, and 
the anisotropic EoS of the magnetized matter. No doubt 
this can be a non trivial, but very interesting and impor- 
tant task. Even though we do not agree with the criticism 
of our paper raised by Comment [l| for the reasons ex- 
plained above, we believe that this discussion may help 
to call attention to an important question about strongly 
magnetized stars that still remains unsolved and needs to 
be investigated self-consistently and in full detail. 
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